Radiating Kerr-Newman black hole in f(R) gravity 
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We derive an exact radiating (nonstationary) Kerr-Newman like black hole solutions to, con- 
stant curvature R = Ro imposed, metric f(R) gravity. This generates a geometry which precisely 
that of radiating Kerr-Newman-de Sitter / anti-de Sitter with f(R) gravity term Ro contributing 
a cosmological-like term. The structure of three horizon-like surfaces, viz. timelike limit, apparent 
horizons and event horizons, are determined. We demonstrate the existence of additional cosmo- 
logical horizons, in f(R) gravity model, apart from regular black hole horizons that exist in the 
analogous general relativity case. In particular, the known stationary Kerr-Newman black hole 
solutions, of f(R) gravity and general relativity, are also retrieved. 
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I. INTRODUCTION 



The f(R) gravity, where f(R) is an analytic function 
of the Ricci scalar R, comes into existence as a straight- 
forward extension of general relativity (GR) Q-Q- I n 
these theories the curvature scalar R of the Lagrangian 
in the Einstein-Hilbert action is replaced by an arbitrary 
function of the curvature scalar thereby modifying GR. 
However, the f(R) action is sufficiently general to en- 
capsulate some of the basic characteristics of higher or- 
der gravity. It is an interesting and relatively simple al- 
ternative to GR, from the study of which yields useful 
conclusions pHl] , including the present day acceleration. 
Unlike GR, which demands metric function derivatives 
no higher than second order, the metric f(R) gravity has 
up to fourth order derivatives This causes compli- 
cations in the calculation and hence, in general, finding 
exact solutions in this theory is laborious. Accordingly, 
little is known about f(R) gravity exact solutions, which 
deserve to be understood better. Nevertheless, recently, 
interesting measures have been taken to get the spher- 
ically solutions of f(R) gravity [5rJl2l]. In particular, 
spherically symmetric black hole (BH) solutions were ob- 
tained for a positive constant curvature scalar in Q , and 
a BH solution was obtained in f(R) gravities by requir- 
ing the presence of a negative constant curvature scalar 
Also, several spherical f(R) BH solutions have been 
obtained 04H1- The generalization of these stationary 
f(R) BHs to the axially symmetric case, Kerr-Newman 
BH, was addressed recently [l3 - [l6l |. In particular, it 



is demonstrated [1J] that the rotating BH solutions for 



f(R) gravity can be derived starting from exact spheri- 
cally symmetric solutions by a complex coordinate trans- 
formation previously developed by Newman and Janis 
fl7j in GR. However, the axially symmetric f(R) is still 



unexplored, e.g., the radiating generalization of the f(R) 
Kerr-Newman BH is still unknown. It is the purpose of 
this paper to generate this metrics and we also present 
the D-dimensional Kerr metric in f(R) gravity. Thus we 
extend a recent work of us [l£| on radiating f(R) BH to 
include rotation. The Kerr metric [l9j is undoubtedly 
the single most significant exact solution in the Einstein 
theory of GR, which represents the prototypical BH that 
can arise from gravitational collapse. The Kerr-Newman 
spacetime is associated with the exterior geometry of a 
rotating massive and charged BH j2pj . 

In this paper, we obtain an exact radiating (nonsta- 
tionary) Kerr-Newman-like metric to, constant scalar 
curvature imposed, f{R) gravity in Section III. We in- 
vestigate further the structure and locations of horizons 
of the radiating f(R) Kerr-Newman metric in Section IV. 
We consider whether f(R) gravity plays any special role 
in the formation of horizons in Section V. The general 
remarks on f{R) gravity is given in Section II and we 
discuss briefly basic equations of HD-/(i?) gravity. The 
paper ends with concluding remarks in Section VII. We 
also give the HD Kerr metric, for constant curvature im- 
posed, f(R) gravity, in the appendix. 

We use units which fix the speed of light and the grav- 
itational constant via G = c = 1, and use the metric 
signature (+, -, -, -). 



II. BASIC EQUATIONS OF METRIC f(R) 
GRAVITY 



In this section we briefly review the metric f(R) grav- 
ity in higher-dimensional (HD) spacetime. The starting 
point is the modified Einstein-Hilbert, D-dimensional, 
gravitational action Q: 



/ = 



f(R)'h 



(i) 
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where g is the determinant of the metric g a b, (a, b — 
Q, 1, D — 1), i? is the scalar curvature, and f(R) is the 



2 



real function defining the theory under consideration. As 
the simplest example, the Einstein-Hilbert action with 
cosmological constant Af is given by f(R) = — (-D— 2)A/. 

From the above action, the equations of motion in the 
metric formalism are just [3j: 



1 



R ab (l + f'(R))--(R + f(R))g a b 
{g a iN 2 - V„V 6 )/'CR) = 2T Qb , 



(2) 



where R a b is the usual Ricci tensor, the prime in f'(R) de- 
notes differentiation with respect to R, and V 2 = V a V a 
with V the usual covariant derivative. Here, we are inter- 
ested in obtaining the constant scalar curvature solutions 
R = Rq. Taking the trace in the Eq. $2$, we get 



2(1 + f'(Ro)) Ro-D (R + f(Ro j) = 0, 



(3) 



where we assume that T = T£ = and note that 
9 a = $ a = D. Eq. determines the negative con- 

stant curvature scalar as [7j 



Ro 



Df(Ro) 



2(1 + f'{Ro)) - D 



DA, 



(4) 



Thus any constant curvature solution R — Rq with 1 
f'(Ro) ^ obeys 1 



Rab — A/gab 



1 + f'(Ro) 



T a b 



(5) 



For this kind of solution an effective cosmological con- 
stant may be defined as A/ = Rq/D for £>-dimensional 
spacetime. In this paper we consider the case with con- 
formal matter (T = T° = 0). For the case of conformal 
matter with non- vanishing Af we have again constant 
R = i?o with R = DAf and g a b is a solution of f(R) pro- 
vided that once again f{DA D ) = A D (2 - D + 2f'DA D )). 
In this case the solution is dS or (A)dS depending on the 
sign of Rq, just as in GR with a cosmological constant. 

We note that the condition 1 + f'(Ro) > required 
to avoid the appearance of ghosts [l( and a necessary 
condition that the stationary f(R) BH becomes a type 
of Schwarzschild-(A)dS BH. Further, we require that 
f"(R) > to avoid the negative mass squared of a scalar- 
field degree of freedom, i.e., to avoid a tachyonic instabil- 
ity [E El- Also f(R) is a monotonic increasing function 
with -1 < f'(R) < 0. 



III. RADIATING KERR-NEWMAN BLACK 
HOLE IN f(R) THEORIES 

Let us begin with the action for f(R) gravity with a 
Maxwell term in the 4D case fill: 



1 

16^ 



-<J 



R + f(R) - F ab F 



ab 



(0) 



The Maxwell tensor is F ao = d a Ab—dbA a , where A a is the 
vector potential. From the variation of the above action 



(|6l), the Einstein equation of motion can be written as 
7, 111 



ab{l + f'(R)) - l(R + f(R)) 9 ab 

9abV 2 -V a V b )f\R) = 2T^ b M , 



with the EMT for charged null dust [2l| 



9ab 



Tab 1 = C(«, r)n a n b + F ap F b " - ^F pa FP° 



(7) 



(8) 



Once again, it is easy to see that the trace of the EMT 
is T EM = 0, due to the fact that F£ = in 4D. How- 
ever, in HD T EM ^ 0. On the other hand, the Maxwell 
equations take the form V a F ab — 0. To proceed further, 
with constant curvature constant Rq, and take the trace 
of the Eq. ([7]) , after some algebra, leads to 



_ 2f(R ) 
R °-7iRoJ—l= 4 ^ 

Thus any constant curvature solution R = Rq obeys 

2 



(9) 



-Rah — Aig ab 



i + /W 



-T, 



ab • 



(10) 



and an effective cosmological constant may be defined as 
A/ =JZo/4 

We are looking for the exterior metric for the gravita- 
tional field produced by the charged radiating and rotat- 
ing object in f(R) gravity, i.e., radiating Kerr-Newman 
BH solutions in f(R) gravity. In retarded time coordi- 
nates the metrics can be expressed as: 



ds 2 



A 



[A - 9a 2 sin 2 6] dv 2 + 2\/~A [dv - a sin 2 d(f\ dr 



.|d02 + A— [A(r 2 + a 2 ) - 8] sin 2 9dvd4> 
4 [A(r 2 + a 2 ) 2 ~ <da 2 sin 2 9} sin 2 ^ 2 



(ii) 



where 




S 2 = 


r 2 + a 2 cos 2 9, 


A = 


r 2 + a 2 -2M(v)r + Q 2 (v) 


e = 


l + ^La 2 cos 2 9, 
12 


A = 




Q 2 (v) = 


2Q 2 (v) 


i + /W 



Apr 2 
12 



(r 2 +a 2 ), 



and the related electromagnetic potential is 



A n 



Q(v)r 
S 2 



VA, 0, 0, -v^asin 2 ^ 



(12) 



Here M(v) and Q{v) are functions of retarded time v 
identified respectively as mass and charge of spacetime, 
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and a is the angular momentum per unit mass. The 
metric (fTTj) is a general solution of the modified Ein- 
stein field equations (fTO]) for a charged null fluid defined 
by the stress energy tensor ijHJ). ft describes the exte- 
rior field of the radiating rotating charged body. Thus 
we have kind of charged radiating rotating metric in de 
Sitter/ anti-de Sitter (dS/AdS) like spacetime or radiat- 
ing Kerr-Newman dS / AdS like solution and we call it 
radiating f(R) Kerr-Newman BH. The stationary f(R) 
Kerr-Newman BH [l5|, 16] in (t, r, 9, </>) can be obtained 
by means of the local coordinate transformations and re- 
placing M(v) and Q{v) by constants M and Q. 

The metric (| 1 1 [) of radiating f(R) Kerr-Newman BH is 
a natural generalization of the stationary Kerr-Newman 
BH solutions of f(R) gravity QJ, [lj], but it is Petrov 
type-II, whereas latter is of Petrov type-D. In addition, 
if Rq = then metric (|1 1JI makes Kerr-Newman metric. 
Hence, we refer solution as radiating f(R) Kerr-Newman 
solution representing gravitational collapse of a charged 
null fluid in a non-flat dS/AdS like spacetime. Thus, the 
metric (| 1 1 [) bears same relation Kerr-Newman as does 
Vaidya metric to Schwarzschild metric. Also for Q = 0, 
the metric (fTTj) is radiating f(R) Kerr spacetime (l5lll6j. 
If in addition i?o — > 0, we have Kerr spacetime and for 
a — > the metric (jTTJ) is f(R) Bonnor- Vaidya spacetime 
which has zero angular momentum (l8| . 

The radiating rotating charged solution discussed here 
are derived under the assumption of the constant curva- 
ture R = Rq, in which case the f(R) models are equiva- 
lent to GR and a cosmological constant and also solution 
is time dependent. 



IV. SINGULARITY AND PHYSICAL 
PARAMETERS OF RADIATING f{R) 
KERR-NEWMAN BH 



invariants are regular everywhere except at the origin 
r = but only at the equatorial plane 9 — tt/2 for 
a, M(v), and Q(v) ^ 0. Hence, the spacetime has the 
scalar polynomial singularity [22j at r — 0. The study of 
causal structure of the spacetime is beyond the scope of 
this paper and will be discussed elsewhere. 

In order to further discuss the physical nature of radi- 
ating f(R) Kerr-Newman BH, we introduce their kine- 
matical parameters. Following [23M28| . the null-tetrad of 
the metric (|11|) is of the form 



In 



VI, 0, 0, -Via sin 2 9 



\A4,— , 1, 0, VA,- a sin 2 , 



V2p 
a 



AiasinO, 0, 



Ai(r 2 + a 2 ) sin 



e' 



Aiasin9, 0, — , VAi(r 2 + a 2 ) sin( 



where 



p = r + ia cos 9 




and p and a are complex conjugates of, respectively, p 
and a. The null tetrad obeys null, orthogonal and metric 
conditions 

l a l a = n a n a = m a m a = 0, l a n a = 1, 
l a m a = n a m a — 0, m a m a = — 1, 

9ab = la n b + kn a - m a fhb - Tnbfh a , 



g- = l a n b + l b n a - m a fh b -m b m a . (14) 



,.nb 



The metric of the radiating f(R) Kerr-Newman BH 
solution has the form pT|) with electromagnetic poten- 
tial given by (|12[) and the energy momentum tensor ([5]l. 
Here, we shall discuss the singularity structure of radi- 
ating f(R) Kerr-Newman BH derived in the previous 
section. The easiest way to detect a singularity in a 
spacetime is to observe the divergence of some invari- 
ants of the Riemann tensor. We approach the singu- 
larity problem by studying the behavior of the Ricci 
R = R a bR ab , (R a b the Ricci tensor) and Kretschmann 
invariants K = R a bcdR abcd , (Rabcd the Riemann tensor). 
For the metric (fTTj) they behave as:- 



R : 



K 



T(Q(v),Ro) 
(r 2 + a 2 cos 2 9) 4 ' 
G{M(v), Q(v), a, cos6>, R Q ) 



(r 2 + a 2 cos 2 



(13) 



where J- and Q are some functions. It is sufficient to 
study the Kretschmann and Ricci scalars for the investi- 
gation of the spacetime curvature singularity (ies). These 



Following Ref. [23j, |24|, a null- vector decomposition of 
the metric is of the form 



g a b = -njb - l a m + 7ab, 



(15) 



where "f a b — m a rhb + mbfh a . The optical behavior of null 
geodesies congruences is governed by the Raychaudhuri 
equation [24l428j . 



de 



1 



— = kq - R ab i a i b - -e 

dv 2 



(16) 



with expansion 0, twist w, shear er, and surface gravity 
K. The expansion [2~i| of the null rays, parameterized by 
v, is given by 



e = vj a 



(17) 



where V is the covariant derivative. In the present case, 
the surface gravity [HJ is 



n a l b V b l a , 



(18) 
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FIG. 1: A plot showing positive roots of Eq. (J2TJ. LEFT: for different values of R and, M(v) = 0.9, Q{v) = 0.3, 9 = 0.3 
and a — 0.3 which correspond to TLSs of radiating f(R) Kerr-Newman BH. It is seen that for Ro > 0, we have three positive 
roots which correspond to r^ LS (inner) and r^ LS (outer), and rdTLS (dS TLS). While for Ro = 0, i.e., in GR limit, have just 
two outer and inner TLSs of Kerr-Newman BHs. RIGHT: for different values of a and, M(v) = 0.9, Q(v) = 0.3, 9 = 0.3 and 
Ro = 0.3 which correspond r^ LS (inner) and r^ LS (outer), and t cltls (dS TLS) of radiating f(R) Kerr-Newman BH. 



and the shear takes form 

(T ab = ®ab - ®Khab- (19) 

The luminosity due to loss of mass is given by Lm — 
—dM/dv, Lm < 1 , and due to gauge charge by Lq — 
—dQ/dv, where Lm,Lq < 1. Both are measured in the 
region where d/dv is timelike [24T - [26l | . 

V. THREE KINDS OF HORIZONS 

A BH has three horizon- like surfaces 24]: timelike 
limit surface (TLS), apparent horizons (AH) and event 
horizons (EH). For a classical Schwarzschild BH (which 
does not radiate), the three surfaces EH, AH, and TLS 
are all identical. Upon switching on the Hawking evapo- 
ration this degeneracy is partially lifted even if the spher- 
ical symmetry stays, e.g., for Vaidya radiating BH, we 
have then AH=TLS, but the EH is different If we break 
spherical symmetry preserving stationarity (e.g., Kerr 
BH), then AH=EH but EH ^ TLS. In general, e.g., for 
radiating Kerr-Newman BH, the three surfaces AH ^ 
TLS ^ EH and they are sensitive to small perturbations. 

Here we are interested in these horizons for the radi- 
ating f(R) Kerr-Newman BH. As demonstrated first by 
York [2J|, the horizons may be obtained to 0(Lm,Lq) 
by noting that (i) for a BH with small dimensionless ac- 
cretion, we can define TLS's (quasi-static limit surface) 
as locus where g{d v , d v ) = g vv = (ii) apparent horizons 
are defined as surface such that 9 ~ and (iii) event 
horizons are surfaces such that dQ/dv ~ 0. 



As mentioned above, in the limit a — > 0, one gets f(R) 
Bonnor- Vaidya solution [l8| , and Eq. (f2"Tj) takes the form 

— r 4 -r 2 + 2M(v)r~Q 2 (v) = 0. (22) 



A. Quasi-static limit surface or TLS 

First, we calculate location of TLS surface, which for 
the nonstationary radiating f(R) Kerr-Newman metric 
requires that prefactor of the dv 2 term in metric vanish; 
It follows form Eq. (JTT]) that TLS will satisfy [H| 

A r - Aga 2 sin 2 9 = 0. (20) 



This equation can be rewritten as 

f(v,r,9) = ^-(l-^fy^+2M(v)r 

+ ^ cos 2 6 sin 2 <9a 4 - cos 6 2 a 2 -Q 2 (v) 42B) 



Equation (f2~Tj) is a reduced quartic equation. It is easy 
to check, under condition of the discriminant in [29j . 
Eq. (|2ip admits four real roots. For positive curvature 
i?o > 0, Eq. (|2"T1) . subject to restriction [29(, has for all 
four real roots as shown in Figures ^ and @. One of 
roots is negative. The other three positive roots corre- 
sponds to r^ LS (inner) and r^ LS (outer), and rdTLS (dS- 
like TLS). Clearly, r^ LS < ^tls < r dTLS and that 
r TLS an d r TLS are TLSs of a BH, whereas the root t^tls 
is supplementary TLS due to the f(R) gravity term. 

I 

This coincides with the nonrotational case in which case 
the various horizons are identified and analyzed by us in 
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FIG. 2: A plot of the function f(v,r,6) for M(v) = Xv and 
Q(v) — fi 2 v 2 with parameters values A = 0.04, [i — 0.08, 8 = 
0.3, a = 0.7 and R = 0.3 



18} and hence, to conserve space, we shall avoid the repe- 
tition of same. Further, in the GR limit Ro — > 0; Q(v) — > 
Q(v), we obtain 

r 2 + cos 2 (0) a 2 - 2 M (y) r + Q 2 (v) = 0, (23) 

which trivially solves to 

r - LS = M(v) - y/M 2 (v)-a 2 cos 2 9-Q 2 (v), 
r+ LS = M(v) + y/M 2 (v)-a 2 cos 2 9-Q 2 (v). (24) 

These are regular outer and inner TLSs for a radiating 
Kerr-Newman BH j2fj, and further in the non- rotating 
limit a — > 0, the solutions (1241) reduces to 



' TLS 



TLS 



= M(v) - ^M 2 {v)-Q 2 (v), 
= M(v) + y / M 2 (v)~Q 2 (v), 



which are TLS of Bonnor-Vaidya BH. Thus the radiating 
f(R) Kerr-Newman BH, in the GR limit and a 0, 
degenerates to Bonnor-Vaidya BH [2l| . 



B. Apparent Horizon 

The AH is the outermost marginally trapped surface 
for the outgoing photons. The AH can be either null or 
space-like, that is, it can 'move' causally or acausally 24 1. 
The AHs arc defined as surfaces such that & ~ 24]. 



Using Eqs. (|14p and (fTS)) give the expression for surface 
gravity 



1 

2S 



dA r 2r 



(26) 



which on inserting the expression for A r , becomes 



Ro 


/ Ro 


12S 2 




J* 


R n a 2 


12 






Eqs. CHJ}, C3 and (J2B) then yield 
r 



12 



R a 2 
12 



+2M(u)r- (a 2 + Q 2 (w) 



(28) 



It is evident that the AHs are zeros of = 0. From 
Eq. (|2"5)l . thus the AH's are given by zeros of 



, s -Ro 4 



12 



(25) 



r 2 +2A-f(w)r-(^a 2 +Q 2 (v) 

m (29) 

There exist, subject to condition 29J, three positive roots 
for Ro > as shown in the Figures ([3]) and (HJ. Unlike, 
TLS, the AH's are 6 independent. Hence, unlike non- 
rotating BHs, they do not coincide in the rotating case. 
The three roots correspond to inner and outer BH AHs, 
and dS-like AH. 



Again in GR limit, we get 



r 2 - 2M(v)r + (a 2 + Q 2 (v)) = 0, 



which admit solutions 



AH 



M(v) - ^M 2 (v) - a 2 - Q 2 (y), 



r+ H = M(y) + y/M 2 (v)-a 2 -Q 2 (y). (31) 

These are regular outer and inner AHs for a radiating 
(30) Kerr-Newman BH, and further in the non-rotating limit 
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FIG. 3: A plot showing positive roots of Eq. ([29]). LEFT: for different values of R and, M(v) = 0.9, Q{v) = 0.3, 6 = 0.3 
and a — 0.3 which correspond to TLSs of radiating f(R) Kerr-Newman BH. It is seen that for Ro > 0, we have three positive 
roots which correspond to r^ H (inner) and r\ H (outer), and t^ah- While for Ro = 0, i.e., in GR limit, have just two outer 
and inner AHs of Kerr-Newman BHs. RIGHT: for different values of a and, M(v) = 0.9, Q(v) = 0.3, 9 = 0.3 and Ro = 0.3 
which corresponds r^ H (inner) and r\ H (outer), and TdAH (dS AH) of radiating f(R) Kerr-Newman BH. 



a 0, the solutions (|31[) correspond to Bonnor-Vaidya 
AHs. Further, Eq. (|3Tj) in the limit a — > becomes ex- 
actly Eq. ([23]). Thus AHs coincide with TLSs for the 
nonrotating but radiating, Bonnor-Vaidya case (18| . 

The discussion in above two subsections are also valid 
for the stationary case discussed in Ref. [Hj]. In the 
stationary case M and Q are constant whereas in the 
radiating case M(v) and Q(v) are function of the retarded 
time v. Thus, Eqs. ([2"Tj) and ([2"5]) are same as derived for 
the corresponding stationary case when M(v) = M 
and Q(v) — Q with M and Q constants. 



C. Event Horizon 

The EH is a null three-surface which is the locus of out- 
going future-directed null geodesic rays that never man- 
age to reach arbitrarily large distances from the BH and 
behave such that d8/dv ~ 0. They are determined via 
the Raychaudhuri Eq. (fl~6]) to 0(Lm,Lq). This defini- 
tion of the EH requires knowledge of the entire future of 
the BH. Hence, it's difficult EH exactly in non-stationary 
spacetime. However, York [24| . gave a working definition 
of the EH, which is in 0(Lm, Lq) equivalent to that the 
photons at EH unaccelerated in the sense that 



dn 2 



~~T~o =0, 



(32) 



with d/dn = n a V a . This criterion enables us to distin- 
guish the AH and the EH to necessary accuracy. It is 
known that (H: 




FIG. 

Q(v) 



A plot of the function g(v,r) for M(v) = Xv and 

0.04, fi = 0.08, 6 = 



fj, 2 v 2 with parameters values A 



0.3 and R = 0.5 



1 



dn 2 



A2Y 2 



l 1 i 2\®^r A r 



2S 



(33) 



For low luminosity, the surface gravity k can be evaluated 
at AH and the expression for the EH can be obtained to 



0{L M ,L Q ). Eqs. 



(|2"rj)) . and the expression for A r , 
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lead to 

Ro 4 ( \ Ro a ' 

— r — 1 

12 V 12 



r 2 + 2A/*(v)r- (a 2 +Q* 2 (v)) = 0, 

(34) 



where 



M*(v) = M(v) + 



VAkE 



L 



M 



(r 2 4- n 2 ) 
Q*(v) = Q(v)+ { - 7 f J l q . 

Eq. is the master equation for deciding the EHs of 
radiating f(R) Kerr- Newman BH. It is interesting to note 
the mathematical similarity with its counterpart Eq. (|29[) 
for AHs. However, unlike the AHs, EHs has 9 depen- 
dence as M*{v) and Q*(v) involve E. For statio nary 
BH, Lm — Lq — 0. Thus, unlike the stationary case |15j . 
where AH=EH ^ TLS, we have shown that for radiating 
f(R) Kerr-Newman BH AH ^ EH ^ TLS. Thus the ex- 
pression of the EH is exactly the same as its counterpart 
AH given by Eq. (|29|) with the mass and charge replaced 
by the effective mass M*(v) and charge Q*{v) pa. l28j. 
The region between the AH and the EH is defined as a 
quantum ergosphere [2~ij ]. The GR limit will lead to the 
same expression as (|3Tj) . with M*(v) and Q*(v) instead 
of, respectively, M(v) and Q{v). 



VI. CONCLUDING REMARKS 

In this paper, we have obtained a class of f(R) Kerr- 
Newman BH in constant curvature R = Rq imposed 
metric f(R) gravity, i.e., a radiating f(R) Kerr-Newman 
metric, describing the gravitational field exterior to radi- 
ating, rotating, charged bodies. The f(R) gravity theo- 
ries are designed to produce a time- varying effective cos- 
mological constant, the BH and spherically symmetric 
solutions of interest are likely to represent central ob- 
jects embedded in cosmological backgrounds. It is evi- 
dent from analysis that the f(R) gravity contributes to 
a cosmological-like term in the solutions and they are 
asymptotically dS/AdS according to the sign of Rq, and 
has geometry of the Kerr-Newman dS/AdS. 

The three kinds of the horizon-like surfaces of the ra- 
diating f(R) Kerr-Newman BHs: TLSs, AHs, and EH 
are studied by the method developed by York [24[ to 
0{Lm,Lq) by a null-vector decomposition of the met- 
ric. It turns out that for each of TLS, AH and EH, there 
exist three surfaces corresponding to the three positive 
roots r~ , r + and r<iAH- As before r~ and r + can be 
viewed, respectively, as inner and outer BH horizons, and 
fdAH as cosmological or dS-like horizon. The fourth root 
r , which is negative also corresponds to cosmologi- 
cal horizon [20[ ■ The analysis presented is applicable to 
stationary f(R) Kerr-Newman BHs as well, but AHs co- 
incide with EHs because stationary BH do not accrete, 
i.e., Lm — Lq — 0. However, the three surfaces no more 
coincide with each other in radiating f(R) Kerr-Newman 



BHs. Thus, we have shown that the presence of the f(R) 
gravity term Rq produces a drastic change in the struc- 
ture of these three horizons. Such a change could have 
a significant effect in the dynamical evolution of these 
horizons. 

The relation between GR and any modified theory of 
gravity is very good way to know how much new theory is 
different from GR. Obviously, When f(R) = R, the the- 
ory reduces to GR. Since for the energy momentum ten- 
sor ©, the trace T = 0, consequently R, f(R) and, f'(R) 
are constant and the theory is equivalent to GR with a 
cosmological constant A/ = Rq/4. Also, the metric f(R) 
gravity corresponds to Brans-Dicke (BD) theory with the 
potential term, Vf» = / - RfiR), = 1 + f'(R) and 
the BD parameter ujb d =0 [l|, |4| . The prototype of gen- 
eral f(R) gravity (non-constant curvature) is the model 
f(R) — R— H 4 /R [J 3) where \i is a mass scale of the 
order of the present value of the Hubble parameter. In 
this model, for large values of Ricci scalar curvature (for 
BH solutions), the f(R) function tends to f(R) = R, and 
kicks in only late as R — > 0; thus the gravity is not linear 
and modified. Thus it admits the solution derived above, 
but with the A/ modified for this choice of f(R). 

To conclude, It is notable that there is no exact solu- 
tion in f(R) gravity coupled to matter with the exception 
of Maxwell fields [11( . We have obtained an exact radiat- 
ing rotating BH in, constant curvature, f(R) gravity for 
charged null dust matter. The solutions presented here 
provides a necessary grounds to further study geomet- 
rical properties, causal structures and thermodynamics 
of these BH solutions, which will be subject of a future 
project. Further generalization of such solutions in more 
general f(R) gravity theories is an important direction 
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Appendix A: Kerr Black holes in HD-/(i?) gravity 

Lastly, we consider a Z?-dimensional Kerr solution in 
constant curvature f(R) gravity, representing generaliza- 
tion of the exterior metric 0] for the axially symmetric 
rotating objects in f(R) gravity. However, we shall re- 
strict ourselves to uncharged case only because if D > 4, 
the trace of the electromagnetic EMT is not zero, which 
is necessary for finding the constant curvature solutions 
from f(R) gravity. 

The Kerr metric describing the stationary, axis sym- 
metric Z?-dimensional spacetime in f(R) gravity may by 
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cast in the form 



ds 2 



A - 9a 2 sin 2 OA dt 



E _ d r 2 - SrdOl 



1 



"6(r 2 



,2\2 



2a 



[0(r 



a e 1 s 

A] sin 2 Oidtdel - r 2 cos 6»id6» 2 - r 2 cos 6>i sin 2 d6> 



A sin 6»i sin 2 B x dSl 



-r 2 cos 6*i sin Q\ 



sin O n dd n+1 , 



(Al) 



where 
E = r 2 + a 2 cos 2 6»i 
1 - 



A = (r 2 +a 2 ) 

e = i 



2R r 2 



(n + l)(n + 2)(n + 3) 
2i? a 2 



(n+l)(» + 2)(n + 3) 



cos t/i. 



D dimensional Kerr-like solutions in dS/AdS spacetime 
and we call it HD-/(i?) Kerr solution. For D = 4, the 
metric ([All) reduces to the Kerr metric in f(R) grav- 
2M ity [UlU In addition, if R = then metric (|A"Tj) 

%r n-2 ' turns out to be the 4D Kerr metric. The correspond- 
ing D-dimensional radiating Kerr metric in f(R) grav- 
ity can be obtained by local coordinate transformations 
(t, r, 9, 4>) — > (v, r, 9, <f) via 



The _D-dimensional solution (|Al|) is a general solution of 
the modified Einstein field Eq. ([2]). The solutions are 
similar to D-dimensional Kerr-dS/AdS solutions. Hence, 
we conclude that the above rotating D-dimensional so- 
lutions of the constant curvature f(R) gravity, is just 



r 2 + a 2 , 
dv — > dv H — — dr, 

dd> — > dd> + —dr. 



(A2) 
(A3) 
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